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Abstract 

The matrix equations of the relativistic random-phase approximation (RRPA) are derived for 
an effective Lagrangian characterized by density-dependent meson-nucleon vertex functions. The 
explicit density dependence of the meson-nucleon couplings introduces rearrangement terms in the 
residual two-body interaction, that are essential for a quantitative description of excited states. Il- 
lustrative calculations of the isoscalar monopole, isovector dipole and isoscalar quadrupole response 
Of208pb, are performed in the fully self-consistent RRPA framework, based on effective interactions 
with a phenomenological density dependence adjusted to nuclear matter and ground-state proper- 
ties of spherical nuclei. The comparison of the RRPA results on multipole giant resonances with 
experimental data constrains the parameters that characterize the isoscalar and isovector channel 
of the density-dependent effective interactions. 

PACS numbers: 21.60.-ii, 21.30.Fe, 21.65.+f, 21.10.-k 



I. INTRODUCTION 



The success of models based on the Relativistic Mean Field (RMF) |T[ approximation 
in describing structure phenomena, not only in nuclei along the valley of /5-stability, but 
also in exotic nuclei with extreme isospin values and close to the particle drip lines, has 
also renewed the interest in theoretical studies based on the Relativistic Random Phase 
Approximation (RRPA) . Although several RRPA implementations have been available since 
the eighties, only very recently RRPA-based calculations have reached a level on which a 
quantitative comparison with experimental data became possible. Two points are essential 
for the successful application of the RRPA in the description of dynamical properties of finite 
nuclei: (i) the use of effective Lagrangians with nonlinear self-interaction terms, and (ii) the 
fully consistent treatment of the Dirac sea of negative energy states. Many studies over the 
last decade have shown that the inclusion of nonlinear meson terms in meson-exchange RMF 
models, or nonlinear nucleon self-interaction terms in relativistic point-coupling models, is 
absolutely necessary in order to reproduce ground-state properties of spherical and deformed 
nuclei on a quantitative level. Techniques which enable the inclusion of nonlinear meson 
interaction terms in the RRPA framework, however, have been developed only recently in 
the calculation of the relativistic linear response 0], and in the solution of the RRPA-matrix 
equation [Q. For a quantitative description of excited states, the RRPA configuration space 
must include not only the usual particle-hole states, but also pair-configurations formed from 
occupied states in the Fermi sea and empty negative-energy states in the Dirac sea. Even 
though it was known for a long time that the inclusion of configurations built from occupied 
positive-energy states and empty negative-energy states is essential for current conservation 
and the decoupling of spurious states [Q, only recently it has been shown that the fully 
consistent inclusion of the Dirac sea of negative energy states in the RRPA is essential for a 
quantitative comparison with the experimental excitation energies of giant resonances 0, ^ . 

The RRPA with nonlinear meson interaction terms, and with a configuration space that 
includes the Dirac sea of negative-energy state, has been very successfully employed in 
studies of nuclear compressional modes p, |^, , of multipole giant resonances and of low- 
lying collective states in spherical nuclei , of the evolution of the low- lying isovector dipole 
response in nuclei with a large neutron excess |10|, and of toroidal dipole resonances |]Tl] . 



An interesting alternative to the highly successful RMF models with nonlinear self- 
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interaction terms, is an effective hadron field theory with medium dependent meson-nucleon 
vertices. Such an approach retains the basic structure of the relativistic mean-field frame- 
work, but could be more directly related to the underlying microscopic description of nu- 
clear interactions. In particular, the density dependent relativistic hadron field (DDRH) 
model Ijr^l has been successfully applied in the calculation of nuclear matter and ground-state 
properties of spherical nuclei |T^, and extended to hypernuclei [0], neutron star matter []T3 



and asymmetric nuclear matter and exotic nuclei |T^. Very recently, in Ref. [|T7l we have 
extended the relativistic Hartree-Bogoliubov (RHB) model ]18| to include density dependent 
meson-nucleon couplings. The effective Lagrangian is characterized by a phenomenological 
density dependence of the a, uj and p meson-nucleon vertex functions, adjusted to properties 
of nuclear matter and finite nuclei. It has been shown that, in comparison with standard 
RMF effective interactions with nonlinear meson-exchange terms, the density-dependent 
meson-nucleon couplings significantly improve the description of symmetric and asymmetric 
nuclear matter, and of isovector ground-state properties oi N ^ Z nuclei. This is, of course, 
very important for the extension of RMF-based models to exotic nuclei far from /3-stability 
(description of the neutron skin, the neutron halo, pygmy isovector dipole resonances), and 
for applications in the field of nuclear astrophysics. 

In this work we derive the RRPA with density-dependent meson-nucleon couplings. Just 
as in the static case the single-nucleon Dirac equation includes the additional rearrangement 
self-energies that result from the variation of the vertex functionals with respect to the 
nucleon field operators, the explicit density dependence of the meson-nucleon couplings 
introduces rearrangement terms in the residual interaction of the RRPA. The rearrangement 
contribution is essential for a quantitative analysis of excited states in the RRPA framework. 
In Sec. H we present the formalism of the relativistic RPA with density-dependent meson- 
nucleon couplings, and derive the RRPA equations in the small amplitude limit of the time- 
dependent RMF. The results of an illustrative calculation of multipole giant resonances in 



208 



Pb are analyzed in Sec. III. Section IV contains the summary and the conclusions. 
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II. FORMALISM OF THE RELATIVISTIC RANDOM-PHASE APPROXIMA- 
TION WITH DENSITY-DEPENDENT MESON-NUCLEON COUPLINGS 



The standard density dependent relativistic hadron field (DDRH) model [112 1 for nuclear 
matter and finite nuclei is defined by the Lagrangian density 

1 



L = tp {i^ ■ d — m) -ip + -{daY — -m^a'^ 



- (1 ^ T3) 

- Qai^ai) - g^ij-f ■ uip - gpij-f ■ pf^) - e^)'^ ■ A ^ ip . (1) 

Vectors in isospin space are denoted by arrows, and bold-faced symbols indicate vectors in 
ordinary three-dimensional space. The Dirac spinor ip denotes the nucleon with mass m. rrirj, 
m^, and nip are the masses of the cx-meson, the cu-meson, and the p-meson. g^, and gp 
are the corresponding coupling constants for the mesons to the nucleon. e^/47r = 1/137.036. 
The coupling constants and the mass of the a-meson are treated as free parameters, adjusted 
to reproduce nuclear matter properties and ground-state properties of finite nuclei. VL^'^ , R^'^ , 
and F^'^ are the field tensors of the vector fields a;, p, and of the photon: 

fi/^^ = a'^cu^ - d^'uo^' (2) 
^/.v = Qi^p^ _ d^pf^ (3) 

The meson-nucleon couplings g„, g^, and gp are assumed to be vertex functions of Lorentz- 
scalar bilinear forms of the nucleon field operators. In most applications of the density- 
dependent hadron field theory these couplings are chosen as functions of the vector density 



pv = yjfj.j'^, with = ip'^fj,ip. Alternatively, the couplings could be functionals of the scalar 
density ps = ipip- It has been shown, however, that the vector density dependence produces 
better results for finite nuclei |]12|, and provides a more natural relation between the self- 
energies of the density-dependent hadron field theory and the Dirac-Brueckner microscopic 
self-energies [0. In the present work we choose the vector density dependence for the 
meson-nucleon couplings. 

The single-nucleon Dirac equation is derived by the variation of the Lagrangian ([^) with 
respect to ip 

idt^P, = {a[-zV - V(r, t)]+V{r, t) + (3{m + 5(r, t)) + S^(r, t)} . (5) 



The Dirac hamiltonian contains the scalar and vector nucleon self-energies defined by the 
following relations: 

S{r,t)=gM^i'r,t) (6) 



Vf,ir, t) = g^{py)uj^{r, t) + gp{pv)r ■ Pf,{r, t) + 



(7) 



The density dependence of the vertex functions g^, guj, and gp produces the rearrangement 
contribution to the vector self-energy 



dgu; J u I dgp - dg„ 



ijjijja . 



(8) 



opv op^ dpv 

The inclusion of the rearrangement self-energies is essential for the energy-momentum con 



servation and the thermodynamical consistency of the model [12, 13 



In the time-dependent RMF model [0 one usually neglects the retardation eff'ects for 
the meson fields, and the self-energies are determined at each time by the solutions of the 
Klein-Gordon and Poisson equations: 



-A 



m, 



-A + mi 



-^Ap{r,t) 



-gaipv) Ps{r,t) 
9Upv)jp{r,t) , 
9p{pv) Jfi{r,t) , 
ejcpir,t) . 



(9) 



This approximation is justified by the large masses in the meson propagators. Retardation 
effects can be neglected because of the short range of the corresponding meson exchange 
forces. The explicit solutions of Eqs. (|^) read, respectively. 



a(r,t) 
Pfiir,t) 

with the Yukawa propagators 



g„{p,)D^{r,r')psir',t)d'r' 



gUpv)D^{r,r')jp{r',t)d^r' 
gp{Pv)Dp{r,r')jp{r',t)d^r' , 
e J D,{r,r')j,p{r',t)d\' , 



(10) 



DJr.r') 



I Q-m^\r-r'\ 

An \r — r'\ 



(11) 



where (j) denotes the a, a;, and p mesons, and the photon. 



The sources of the fields are the nucleon densities and currents calculated in the no-sea 
approximation 

A 

1=1 

A 



i=l 

A _ 



i=l 
Z 



Jc>^ir,t) = J2^,ir,th^t^ir,t) . (12) 

i=l 

where the summation is over all A occupied states in the Fermi sea, i.e. only occupied 
single-nucleon states with positive energy explicitly contribute to the nucleon self-energies. 
Even though the stationary solutions for the negative-energy states do not contribute to 
the densities in the no-sea approximation, their contribution is implicitly included in the 
time-evolution of the nuclear system [§, 

The Relativistic Random Phase Approximation (RRPA) represents the small amplitude 
limit of the time- dependent relativistic mean-field theory. In the remainder of this section 
we will derive the RRPA equations with density- dependent meson-nucleon couplings from 
the response of the density matrix p{t) to an external field 

F{t) = Fe"*"* + h.c. , (13) 

which oscillates with a small amplitude. In the single-particle space this field is represented 
by the operator 

fit) = E Ut) (14) 

kl 

The expression for the single-particle density matrix reads 

/3(r,r',t)=^|^,(r,t))(^,(r',t)| . (15) 

i=l 

By writing the Dirac spinor in terms of large and small components 

Mr,t))=( , (16) 

the density matrix takes the form 



p(r,r ,t) 



i=l i=l 
A ^ A 



^^E9^{r,t)f|{r',t) E 9i{r,t)gj{r' ,t) , 

\ i=l 1=1 / 
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(17) 



The equation of motion for the density operator reads 



idtp = h{p) + f{t),p 



and in the small amplitude limit the density matrix is expanded to linear order 



Pit) = p(°) + 6p{t) 



(18) 



(19) 



where p^^^ is the stationary ground-state density. From the definition of the density matrix 
(|T5|), it follows that p{t) is a projector at all times, i.e. p{t)^ = pif)- In particular, this 
means that the eigenvalues of p'-^^ are and 1. In the non-relativistic case particle states 
above the Fermi level correspond to the eigenvalue 0, and hole states in the Fermi sea 
correspond to the eigenvalue 1. In the relativistic case, one also has to take into account 
states from the Dirac sea. In the no-sea approximation these states are not occupied, i.e. 
they correspond to the eigenvalue of the density matrix. 



(0) 

Pki 



(0) 



for unoccupied states above the Fermi level (index p) 

1 for occupied states in the Fermi sea (index h) (20) 
for unoccupied states in the Dirac sea (index a) . 



p{tf' = p (t) also implies, in leading order. 



p(°)5p + 5pp(°) = 6p 



(21) 



and this means that the only non-vanishing matrix elements of 6p are: Spph, Sphp, Sp^h, and 
6pha- These are determined by the solution of equation (|TH]), which in the linear approxi- 
mation reads 



idtSp= \h^^\6p 



+ 



dh 
dp 



5p,p 



(0) 



+ 



(0) 



with 



dh dh dh ^ ^ dh dh 

^P ph ^Pph ^Php ah '^P'^h OPha 



(22) 



(23) 



Under the influence of the external field (p!3|), in the small amplitude limit 5p also exhibits 
the harmonic time dependence e"*"^*. Taking into account that h^^^ = dki^k is diagonal in 
the stationary basis, the resulting RRPA equations read 



{u — ep + eh)Spph — fph + X! ^ph'hp'Spp'h' + Vpp'hh'Sph'p' + X! ^ph'ha'SPa'h' + Vpa'hh'^Ph', 

p'h' a'h' 
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{uj - ea + eh)Spah 
{uj -eh + ep)6php 

{uj - eh + ta)5pha 

or, in matrix form 



fah + X! ^ah'hp'Spp'W + Vap' hh'S Ph'p' + ^ Vah'ha'^ Pa' h' + Vaa' hh'^ Ph' a' 
p'h' a'h' 

fhp + X! ^hh'pp'SPp'h' + Vhp'ph'Sph'p' + X! ^hh'pa'SPa'h' + Vha'ph'^ Ph'a' 
p'h' a'h' 

fha + X] ^hh'ap'Spp'h' + Vhp'ah'S Ph'p' + X! ^/i/!.'aa"5Pa'h' + Vfta'aft,' V/i'a' 

(24) 



1 
-1 



A 

B* A* j 



X 



F 
F 



(25) 



The RRPA matrices A and B read 



A 



B 



^h)^pp'^hh 



Vpp'hh' Vpa'hh' 
y Vap'hh' Vaa'hh' 



{^a — ^h)8aa'5. 



hh' 



Vph'hp' Vph' ha' 
Vah'hp' Vah' ha' 



and the amphtudes X and Y are defined 

/ 

X = 



5pph 

\ 5 Pah 



Y 



5php 
5pha 



The vectors which represent the external field contain the matrix elements 





fph 




fhp \ 


F = 










y fah y 




y fha J 



(26) 
(27) 



(28) 



(29) 



In the self-consistent RRPA the matrix elements of the residual interaction are derived from 
the Dirac hamiltonian of Eq. (|^) : 

dhnr 



abed 



dp 



'db 



(n)^6+(r-2)V(ri, r2)^c(n)^d(r'2)rfVic?V2 . 



(30) 



In order to calculate the contributions of each meson field to V{ri,r2), we expand the 
meson-nucleon couplings and their derivatives around the ground-state density p° 

dgi 



9i{pv) = 9i{pl) + 
dpv dpv 



dpv 
' 9pl 



Spv 



5pv 



(31) 
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If for the meson fields appearing in the scalar and vector nucleon self-energies we use the 
explicit solutions ( [T0| ) in terms of the meson propagators and nucleon densities and currents, 
the individual contributions of the meson fields to V{ri, are obtained from the particle- 
hole matrix element of the Dirac hamiltonian: 

• the contribution of the isoscalar-scalar sigma field 

V^{ri,r2) = -f3if329a{pviri))g^{pvir2))D^{ri,r2) 

[ dp^[ri) dpliri) dp^{ri) J ri 

- \(3il2ga{pv{ri)) . Ps{r2) + ti(32 r.^^,'' . ps{ri)ga{pv{r2)) 

'''' 

where 

lairi) = J rg^{py{r))Dl{ri,r)ps{r)dr ; 

• the contribution of the isoscalar-vector omega field 

Vu,{ri,r2) = {Pl^)l{Pl^,)2gUPv{rl))g^{p^{r2))D^{r^,r2) 

i dpv{ri) dpi{ri) J ri 

+ <gu^{pv{ri)) ^ , . Pv{r2) + o / \ P 
[ opv{r2) Opv[ri) 

+ rf^i" . pv{ri) p^(r2)l]lil2^^(r'i,-r2) , (33) 

dp^{ri) dp^{r2) J 

where 

Iwiri) = J rg^{p^{r))D^{ri,r)p^{r)dr ; 

• the contribution of the isovector-vector rho field 

V/Cri.r,) = (/37")i(/J7p)2TMg,(p„(r,))9,(().(''2))B,(n,r2) 



\dpy{ri) ^ dp^{ri) ^ dpl{ri) " j ri 
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where 

Ip{ri) = I rgp{py{r))Dl{ri,r)pty{r)dr 



• and finally the contribution of the Coulomb field 

V,in,r2) = e\(3Y)i{Pl,)2D,{r^,r2) . (35) 

the subscripts 1 and 2 of the Dirac matrices refer to particle 1 and 2, respectively, p^,, ps and 
Ptv denote the vector, scalar, and isovector-vector density, respectively, and the derivatives 
of the meson-nucleon couplings with respect to the vector density are evaluated at ground- 
state density p°. The radial integrals (0 = a, u, p) contain D^Ar, r'), which is the radial 



factor in the first term of the multipole expansion of the Yukawa propagator (p 

1 oo L 
' " L=0 M=-L 

We notice that, in addition to the direct contribution of the meson exchange interac- 



tions (first terms in Eqs. (|32D -(|34D), the explicit density dependence of the meson-nucleon 
couplings introduces a number of rearrangement terms in the residual two-body interaction 
y{fiif2)- These rearrangement terms are essential for fully consistent RRPA calculations. 
Only when their contribution is included in the matrix elements of the residual interaction, 
it becomes possible to reproduce reasonably well the excitation energies of giant multipole 
resonances. Without rearrangement terms, one finds discrepancies of the order of several 
MeV between the experimental excitation energies and the RRPA peak energies, calculated 
with relativistic effective interactions that are adjusted to ground state properties of spher- 
ical nuclei. A very similar effect is observed in RRPA calculations based on effective forces 
with non-linear meson self-interactions, when the contribution of the non-linear terms is not 
included in the matrix elements of the residual interaction [Bl. 



III. ILLUSTRATIVE RRPA CALCULATIONS: GIANT RESONANCES 

In this section the RRPA with density-dependent meson-nucleon couplings is applied in 
illustrative calculations of giant resonances in spherical nuclei. In particular, we analyze the 
isoscalar monopole, the isovector dipole and the isoscalar quadrupole resonances in ^°^Pb. 
We will show which isoscalar and isovector properties of the effective mean-field interactions 
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affect the multipole strength distributions, and how the results of RRPA calculations can 
be used to constrain the effective interaction. 

For the density dependence of the meson-nucleon couplings we adopt the functionals used 
in Refs.[]TB|, |T^, |T^. The coupling of the a-meson and u;- meson to the nucleon field reads 

9i{p) = gt{.Ps^t)fi{,x) for i = a,uj, (37) 

where 



/^(^) = ;^T. 2 (38) 



l + hi{x + diY 
'1 + Ci{x + di 

is a function of x = p/psat, and psat denotes the baryon density at saturation in symmetric 
nuclear matter. The eight real parameters in (^) are not independent. The five constraints 
/j(l) = 1, /"(I) = /"(I), and /j"(0) = 0, reduce the number of independent parameters to 
three. Three additional parameters in the isoscalar channel are: ga{Psa.t), 9u}{Psa.t), and m„ 
- the mass of the phenomenological sigma-meson. For the p-meson coupling the functional 
form of the density dependence is suggested by DB calculations of asymmetric nuclear 
matter 

9p{p) = ^p(Psat) exp [-ap{x - 1)] . (39) 

The isovector channel is parameterized by gp{psa.t) and a^. For the masses of the u and p 
mesons usually the free values are used: = 783 MeV and rup = 763 MeV. In principle one 
could also consider the density dependence of the meson masses. However, since the effective 
meson-nucleon coupling in nuclear matter is determined by the ratio g/m, the choice of a 
phenomenological density dependence of the couplings makes an explicit density dependence 
of the masses redundant. 

The eight independent parameters, seven coupling parameters and the mass of the a- 
meson, are adjusted to reproduce the properties of symmetric and asymmetric nuclear mat- 
ter, binding energies, charge radii and neutron radii of spherical nuclei. In particular, in 
Ref. fl^ we have introduced the density-dependent meson-exchange effective interaction 
(DD-MEl), whose parameters are displayed in Table |l|. The seven coupling parameters 
and the a-meson mass have been simultaneously adjusted to properties of symmetric and 
asymmetric nuclear matter, and to ground-state properties (binding energies, charge radii 
and differences between neutron and proton radii) of twelve spherical nuclei. For the open 
shell nuclei pairing correlations have been treated in the BCS approximation with empirical 
pairing gaps (five-point formula). 
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In Ref. the relativistic Hartree-Bogoliubov (RHB) model with the density-dependent 
interaction DD-MEl in the p/i-channel, and with the finite range Gogny interaction DIS in 
the pp-channel, has been tested in the analysis of the equations of state for symmetric and 
asymmetric nuclear matter, and of ground-state properties of the Sn and Pb isotopic chains. 
It has been shown that, as compared to standard non-linear relativistic mean-field effective 
forces, the interaction DD-MEl has better isovector properties and therefore provides an 
improved description of asymmetric nuclear matter, neutron matter and nuclei far from 
stability. 

In the present analysis we perform fully consistent RRPA calculations of isoscalar 
monopole, isovector dipole and isoscalar quadrupole giant resonances in ^°^Pb. The single- 
particle basis and the particle-hole couplings are obtained from the same effective La- 
grangian, and the configuration space includes both particle-hole pairs, as well as pairs 
formed from hole states and negative-energy states from the Dirac sea. Our starting point 
is the DD-MEl effective force, both in the Dirac hamiltonian (j^), as well as the residual 
interaction. We then proceed to construct families of density-dependent interactions with 
some given characteristic (compressibility, asymmetry energy, etc.), and study the resulting 
properties of giant resonances. 

The isoscalar giant monopole resonance (ISGMR) represents the most simple mode of 
collective excitations in nuclei. In particular, the ISGMR in heavy nuclei is the only source 
of experimental information on the nuclear matter compression modulus Koo- This quantity 
determines basic properties of nuclei, supernovae explosions, neutron stars and heavy-ion 
collisions. The range of values of K^o has been deduced from the measured excitation 
energies of the ISGMR in spherical nuclei. The presently available experimental data set, 
however, does not limit the range of K^o to better than 200 — 300 MeV. The microscopic 
determination of the nuclear matter compressibility is based on the construction of sets 
of effective interactions that differ mostly by their prediction of the excitation energies of 
ISGMR, i.e. by the value of Kao, but otherwise reproduce reasonably well experimental data 



on ground-state nuclear properties [^, ^ . Effective interactions with different values of K^o 
are used to calculate bulk ground-state properties of heavy spherical nuclei in a self-consistent 
mean-field framework, and RPA or time- dependent mean-field calculations are performed for 
the isoscalar monopole excitations. Such a fully consistent calculation of both ground-state 
properties, as well as ISGMR excitation energies, restricts the range of possible values for 
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KoQ. However, since there are also other effects beyond the mean- field level which influence 
the isoscalar monopole resonance (anharmonicities, pairing, coupling between single-nucleon 
and collective motion), it has been argued that, rather than on the systematics over 
the whole periodic table, the determination of the nuclear compressibility should rely more 
on a good measurement and microscopic calculations of GMR in a single heavy nucleus such 
as ^°*Pb. Microscopic calculations have been performed both in the non-relativistic and 
in the relativistic mean-field framework. Modern non-relativistic Hartree-Fock plus RPA 
calculations, using both Skyrme and Gogny effective interactions, indicate that the value of 
Kcxy should be in the range 210-220 MeV gl], H. In particular, in Ref. fH a set of effective 
Gogny forces was generated, which on one hand allow a good description of static properties 
of nuclei, and on the other hand span the range 210 < K^q < 300 MeV. It was shown that the 
RPA calculations reproduce the available experimental data on ISGMR in medium-heavy 
and heavy nuclei only for Koo in the range 210-220 MeV. In Ref. it has been shown 
that even generalized Skyrme forces, with both density- and momentum-dependent terms, 
can only reproduce the measured breathing mode energies for values of Koo in the range 
215 ± 15 MeV. In relativistic mean-field models based on non-linear meson self-interactions 
on the other hand, results of both RRPA and time-dependent calculations suggest that 
empirical GMR energies are best reproduced by an effective force with Koo ~ 250 — 270 
MeV |2^. It has to be emphasized, however, that even though relativistic calculations 
have been performed using non-linear effective interactions with different values of Koo, 
these forces were not constructed specifically with the purpose of determining Koo. Rather, 
standard non-linear effective interactions have been used, which also exhibit other differences 
that could affect the microscopic determination of the nuclear matter compressibility. 

Starting from DD-MEl, in this work we have generated a consistent set of relativistic 
density-dependent effective interactions with 220 < < 280 MeV. The same functional 
form for the density dependence for the meson-nucleon couplings has been used for these 
forces and, except for the value of Koo, their parameters have been adjusted to the same set 
of experimental data on ground-state properties of twelve spherical nuclei [1^. The results of 
fully consistent RRPA calculations with these forces are shown in Fig. |I], where we display 
the calculated excitation energies of ISGMR in ^°^Pb as function of the nuclear matter 
compressibility. The shaded region denotes the range of presently available experimental 
data [p5|. We notice that, in accordance with the results obtained with relativistic effective 
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forces with non-linear meson self-interactions, only the density-dependent interactions with 
Koo ~ 260 — 270 MeV reproduce the experimental value. We have also verified that ISGMR 
excitation energies for lighter nuclei, calculated with these particular interactions, are closest 
to the empirical curve ~ 80 A^^^^^ MeV and that they reproduce the experimental 
excitation energies. For the density-dependent effective interaction with = 270 MeV, 
in Fig. we display the isoscalar monopole strength distribution and transition densities 
in ^''^Pb. The position of the ISGMR peak is ai E = 14.1 MeV, and we plot the proton, 
neutron and total isoscalar transition densities. 

The present analysis, therefore, confirms that there is a pronounced difference between the 
values of the nuclear matter compression modulus predicted by microscopic non-relativistic 
{K^ ^ 210-230 MeV) and relativistic {K^ ^ 250-270 MeV) mean-field plus random phase 
approximation calculations. The origin of this discrepancy is at present not understood, even 
though there are some indications that it might be due, at least in part, to the differences in 
the density dependence of the asymmetry energy predicted by non-relativistic and relativistic 
models . 



The calculated properties of isovector dipole giant resonances (IVGDR) will be predom- 
inantly determined by the isovector channel of the effective interaction. In particular, the 
excitation energies of IVGDR can be directly related to the nuclear matter asymmetry en- 
ergy. The energy per particle of asymmetric nuclear matter can be expanded about the 
equilibrium density psat in a Taylor series in p and a 



where 



E{p, a) = E{p, 0) + ^2(p)«' + 54(p)«' + ■ ■ ■ (40) 



a^^. (41) 
N + Z ^ ' 



E{p, 0) = -a, + (p - Psat)' + ... (42) 



and 



S2{p) = ^4 + — (P - Psat) + (P - Psat) + " " " (43) 

Psat -LoPsat 

The empirical value of the asymmetry energy at saturation density (volume asymmetry) 
'S'2(psat) = a4 = 30 ± 4 MeV. The parameter po defines the linear density dependence of the 
asymmetry energy, and AKq is the correction to the incompressibility. The contribution of 
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the term 6*4 (p)a;^ in ( ^OD is very small in ordinary nuclei and the coefficient is not constrained 
in the mean-field approximation. 

A ground-state nuclear property which is directly determined by the asymmetry energy 
is the difference between the neutron and the proton radii. In a recent study of neutron 
radii in non-relativistic and covariant mean-field models |]2S|, the linear correlation between 
the neutron skin and the symmetry energy has been analyzed. In particular, the analysis 
has shown that there is a very strong linear correlation between the neutron skin thickness 
in ^°^Pb and the individual parameters that determine the symmetry energy 6*2 (p): CI4, po 
and AKq. The empirical value of r„ — rp in ^°^Pb (0.20 ± 0.04 fm from proton scattering 



data ||29[, and 0.19 ± 0.09 fm from the alpha scattering excitation of the isovector giant 
dipole resonance ||3^) places the following constraints on the values of the parameters of 
the symmetry energy: 04 ^ 30 — 34 MeV, 2 Mev/fm^ < Po < 4 Mev/fm^, and —200 MeV 
< Ai^o < -50 MeV. 

Properties of isovector collective modes in finite nuclei should, in principle, provide ad- 
ditional constraints on the isovector channel of the effective interaction. In an analysis of 



Skyrme forces and giant resonances in exotic nuclei [31|, Reinhard noticed a somewhat sur- 
prising property of the IVGDR: while it is true that the excitation energy of this resonance 
is sensitive to the volume asymmetry 04, the resonance energy decreases by increasing the 
asymmetry energy at saturation. This was qualitatively explained by noticing that an in- 
crease in the volume asymmetry is always accompanied by an increase of the slope po, i.e. 
of the linear density dependence of the asymmetry energy. In order to study this effect in 
a more quantitative way, we have generated, starting from DD-MEl, a set of eight den- 
sity dependent effective interactions with 30 MeV < 04 < 37 MeV. The parameters of the 
density-dependent meson- nucleon couplings have been adjusted in such a way that, while 
increasing 04 in units of 1 MeV, the resulting effective interactions still reproduce the same 
set of data on ground-state properties of spherical nuclei, that was used for the original 
interaction DD-MEl |T^. This means that these effective interactions essentially differ only 
in their description of the asymmetry energy curve as function of the baryon density. 

The resulting nuclear matter asymmetry energy curves, and the calculated IVGDR exci- 
tation energies in ^°'^Pb, are displayed in Fig. |^. In the upper left panel we plot the RRPA 
excitation energy of the IVGDR in ^osp]^ function of the volume asymmetry a^. Similar to 
what has been observed in Ref. the resonance energy decreases with increasing 04. The 
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reason for this decrease is shown in the lower left panel, where we plot the corresponding val- 
ues of the slope parameter po, which defines the linear density dependence of the asymmetry 
energy. We notice that, in order to reproduce the bulk properties of spherical nuclei, an in- 
crease of 04 necessitates a non-linear increase of Pq. The resulting asymmetry energy curves 
as functions of the baryon density are shown in the right panel of Fig. |^. The increase of Pq 
with 04 implies a transition from a parabolic to an almost linear density dependence of 5*2 in 
the density region p < 0.2 fm~^. This means, in particular, that the increase of the asymme- 
try energy at saturation point will produce an effective decrease of 5*2 below p ~ 0.1 fm^^. 
But this is, of course, the density region characteristic for the IVGDR. We find, therefore, 
that the excitation energy of the IVGDR decreases with increasing S'2(psat) = 04, because 
this increase implies a decrease of S2 at low densities characteristic for the surface modes. In 
the upper left panel of Fig. ^ we also compare the calculated IVGDR peak energies for ^°®Pb, 
with the experimental value of 13.3 ± 0.1 [0. It appears that the experimental IVGDR 
excitation energy constrains the nuclear matter asymmetry energy at saturation density to 
the interval 34 MeV < 04 < 36 MeV. For the effective interaction with 04 = 35 MeV, in 
Fig. ^ we display the RRPA isovector dipole strength distribution and the corresponding 
proton, neutron, and total isovector transition densities for the peak at 13.3 MeV in ^°^Pb. 

Fig. ^ illustrates what happens when the increase of the nuclear matter asymmetry energy 
at saturation density is not accompanied by an increase of the slope parameter p^. Starting 
with DD-MEl, which has 04 = 33.1 MeV, we have generated a set of effective interactions 
with different values of 04, but now they all have the same slope parameter po (lower left 
panel), i.e. the parameters are not readjusted to reproduce the data set of ground state 
properties of spherical nuclei. Binding energies and radii are only approximately reproduced 
with these effective interactions. The resulting asymmetry energy curves as functions of the 
baryon density are plotted in the right panel. Since po is constant, by increasing 04 the 
asymmetry energy 5*2 increases for all densities. As a result, the IVGDR peak energies of 
208p^ increase linearly with 04 (upper left panel). 

In non-relativistic RPA calculations, the excitation energy of the isoscalar giant 
quadrupole resonance (ISGQR) can be directly related to the nucleon effective mass that 
characterizes a given effective interaction. In the non-relativistic mean-field approximation, 
the total effective mass m* of a nucleon in a nucleus characterizes the energy dependence 
of an effective local potential that is equivalent to the, generally nonlocal and frequency 



16 



dependent, microscopic nuclear potential |Q. m* is a measure of the density of single- 



nucleon states around the Fermi surface and, therefore, it affects the giant resonances. For 
Skyrme interactions, in particular, a linear dependence on m* is found for the RPA exci- 
tation energies of the ISGQR. The larger the effective mass, i.e. the higher the density of 
states around the Fermi surface, the lower is the calculated ISGQR excitation energy. Both 
the calculation of ground-state properties in spherical nuclei, as well as the RPA results for 
ISGQR excitation energies, place the following constraint on the nucleon effective mass for 



Skyrme-type interactions: m*/m = 0.8 ± 0.1 |31 



In the relativistic framework the concept of effective nucleon mass is more complicated. 
In addition to the k — mass (characterizes the momentum dependence of the mass opera- 
tor) and E — mass (characterizes the explicit energy dependence of the mass operator), a 
third effective mass, the " Lorentz mass" appears in the relativistic approach. It results from 



different Lorentz transformation properties of the scalar and vector potentials p^, 35]. How 



ever, what is usually termed "the relativistic effective mass" that characterizes an effective 



interaction, is a fourth mass: the Dirac mass 34 



mjj = m + S{r) , (44) 

where m is the nucleon mass and S{r) is the scalar nucleon self-energy. It has to be em- 
phasized that the Dirac mass should not be identified with the effective mass determined 
from non-relativistic shell and optical model analyses of experimental data, i.e. with the 
"non-relativistic-type effective mass". The Dirac mass is determined, on one hand by the 
binding energy at saturation density in nuclear matter (the effective single-nucleon potential 
is the sum of the attractive scalar and repulsive vector nucleon self-energies), and on the 
other hand by the empirical spin-orbit splittings in finite nuclei (the effective single-nucleon 
spin-orbit potential is proportional to the difference between the scalar and vector self en- 
ergies). This is the reason why, for virtually all mean-field relativistic effective interactions, 
0.55m < m£, < 0.60m. 

In contrast to the non-relativistic effective mass, the Dirac mass cannot be related to 
the ISGQR. On the other hand, we would like to use the isoscalar quadrupole response to 
constrain the isoscalar properties of our density-dependent effective interactions. This can 
be done in the following way. We first notice that a particular ratio of isoscalar parameters 
bi and Cj {i = a, u) in Eq. (|5B|), characterizes the density of single-nucleon states around the 
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Fermi surface. In Fig. ^ we display the average energy gaps between the last occupied and 
the first unoccupied major shells < -Eparticies > — < -Ehoios > in ^°^Pb, where 

and the sums run over occupied (unoccupied) states within a major shell. The average gaps 
for proton and neutron states are plotted as functions of the parameter 6 

SJ^, (46) 

see Eq. p^). Starting from DD-MEl, we have generated a set of five effective interac- 
tions with 0.93 < 5 < 1.01. For each interaction the remaining parameters were readjusted 
to reproduce our standard set of ground-state data for twelve spherical nuclei, as well as 
the nuclear matter equation of state. The average gap between the last occupied and first 
unoccupied major shells, both for proton and neutron states, is approximately linearly pro- 
portional to 5. This parameter, therefore, plays the role of the inverse of the effective 
mass. As functions of 5, in Fig. |^ we plot the corresponding centroid energy of the isoscalar 
quadrupole Hartree response in ^°^Pb (upper panel), and the peak energies of the ISGQR 
obtained by the full RRPA calculation with the five density-dependent interactions (lower 
panel). The calculated ISGQR excitation energies are compared with the experimental value 
of 10.9 ± 0.3 MeV (shaded area). Both the centroids of the Hartree response and the 
ISGQR peak energies are linearly proportional to 5 and the comparison with experimental 
data on ISGQR, therefore, places an additional constraint on the parameters that charac- 
terize the isoscalar channel of the effective interaction. For 5 = 0.93, in Fig. ^ we plot the 
RRPA isoscalar quadrupole strength distribution in ^°'^Pb (left panel), and for the ISGQR 
peak at 11.2 MeV the proton, neutron, and total isoscalar transition densities. The position 
of the calculated peak should be compared with the empirical excitation energy 10.9 ± 0.3 
MeV, and also the Ohu low-lying discrete 2+ state at 4.62 MeV is found in good agreement 
with the experimental value of 4.07 MeV. 

IV. SUMMARY AND CONCLUSIONS 

During the last decade the standard RMF models with nonlinear meson-exchange effec- 
tive interactions have been very successfully applied in the description of a variety of nuclear 
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structure phenomena. In the last couple of years also the relativistic random-phase approxi- 
mation (RRPA), based on effective Lagrangians with nonlinear meson self- interaction terms, 
has been used to investigate properties of low-lying collective states and of giant resonances. 
The use of nonlinear effective interactions, however, presents not only a number of techni- 
cal problems, but also the predictive power of models based on these type of interactions 
appears to be somewhat limited, especially for isovector properties of exotic nuclei far from 
/3-stability. An interesting alternative are models with density-dependent meson-nucleon 
vertex functions. Even though these two classes of models are essentially based on the 
same microscopic structure, i.e. on density dependent interactions, the latter can be more 
directly related to the underlying microscopic nuclear interactions. In a number of recent 
analyses it has been also shown that relativistic effective interactions, with explicit density 
dependence of the meson-nucleon couplings, provide an improved description of asymmetric 
nuclear matter, neutron matter and nuclei far from stability. 

Among the new structure phenomena observed or predicted in nuclei far from stability, 
one of the most interesting is the evolution of the isovector dipole response in nuclei with 
a large neutron excess. The multipole response of nuclei with large neutron excess has 
been the subject of many theoretical studies in recent years, and some predictions have 
been confirmed by very recent experimental data on low-lying electric dipole strength in 
neutron rich nuclei [ P7| , There are, however, many unknowns and this topic presents 
an interesting challenge for modern theoretical advances. It is, therefore, important to 
develop also a relativistic framework, based on effective Lagrangians with density-dependent 
meson-nucleon couplings, in which the dynamics of exotic collective modes in nuclei far from 
stability can be investigated. 

In this work we have derived the RRPA matrix equations in the small amplitude limit 
of the time-dependent relativistic mean-field theory. The explicit density dependence of the 
meson-nucleon vertices introduces a number of rearrangement terms in the residual two-body 
interaction. We have found that the rearrangement contribution to the matrix elements of 
the RRPA equations is crucial for a quantitative comparison with experimental data on 
giant resonances. In the present analysis we have performed illustrative RRPA calculations 
of the isoscalar monopole, isovector dipole and isoscalar quadrupole response of ^°^Pb. The 
calculations are fully self-consistent: the single-particle basis and the particle-hole couplings 
are generated from the same effective Lagrangian, and the RRPA configuration space in- 
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eludes both the positive-energy particle-hole pairs, as well as pairs formed from hole states 
and negative-energy states in the Dirac sea. On one hand, we have tested our approach by 
comparing the RRPA results for giant resonances with well known experimental data. On 
the other hand we have also analyzed how the RRPA results on multipole giant resonances 
can be used to constrain the parameters that characterize the isoscalar and isovector chan- 
nel of the density-dependent effective Lagrangians. Starting with the recently introduced 
effective interaction DD-MEl Jl^, RRPA calculations have been performed for families of 
density-dependent interactions with a given characteristic (nuclear matter incompressibility, 
asymmetry energy etc.). 

The analysis of the isoscalar monopole response has shown that only the density- 
dependent interactions with the nuclear matter compression modulus in the range Kao ~ 
260 — 270 MeV, reproduce the experimental excitation energy of the isoscalar giant monopole 
resonance in ^°'^Pb. This confirms our previous results obtained with relativistic effective 
forces with non-linear meson self-interactions and points, once again, to the pronounced 
difference between the values of the nuclear matter compression modulus predicted by mi- 
croscopic non-relativistic and relativistic mean-field plus RPA calculations. The RRPA 
results for the isovector dipole response constrain the isovector channel of the effective inter- 
actions. By using interactions with different values of the volume asymmetry energy 04, but 
which otherwise reproduce the same data set of ground-state properties of spherical nuclei, 
we have shown that the calculated IVGDR peak energy actually decreases by increasing the 
asymmetry energy at saturation. The comparison with the experimental IVGDR excitation 
energy constrains the volume asymmetry to the interval 34 MeV < 04 < 36 MeV. In the 
non-relativistic framework the isoscalar quadrupole response can be related to the effective 
mass of the mean-field interaction. The concept of effective mass in the relativistic mean- 
field models is more complicated, and the quantity which is usually termed as "effective 
mass" cannot be identified with the effective mass determined from non-relativistic shell 
and optical model analyses of experimental data. Nevertheless, we have shown that a com- 
parison of RRPA results with the empirical ISGQR and with the low-lying Qhuj 2+ state, 
places an additional constraint on the parameters which characterize the isoscalar channel 
of the density-dependent effective interactions. 

The RRPA with density-dependent meson-nucleon couplings presents an important step 
in the relativistic description of the nuclear many-body problem. In the present work we did 
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not attempt an analysis of the multipole response in exotic nuclei far from /3-stability. In 
order to do that, pairing correlations must be included in the RRPA framework. Work is in 
progress on the fully self-consistent relativistic quasiparticle random-phase approximation 

(RQRPA), based on effective Lagrangians with density-dependent meson-nucleon couplings, 
and formulated in the relativistic Hartree-Bogoliubov canonical single-particle basis. 
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TABLE I: The effective interaction DD-MEl. See Eqs. (|3|)- 
parameters. 
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FIG. 1: Density-dependent RRPA peak energies of the ISGMR in 208p|-, function of the nuclear 
matter compressibiHty K^o- The calculated peaks are shown in comparison with the experimental 
excitation energy of the monopole resonance: E = 14.1 it 0.3 MeV [25]. 



FIG. 2: The isoscalar monopole strength distribution (left panel) and transition densities (right 
panel) in ^'^^Pb, calculated with a density-dependent effective interaction with K^o = 270 MeV. 
The proton, neutron and total isoscalar transition densities correspond to the peak at E = 14.1 
MeV. 
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FIG. 3: The IVGDR excitation energy of ^'^^Pb (upper left panel), and the parameter po of the 
linear density dependence of the nuclear matter asymmetry energy, as functions of the volume 
asymmetry 04. The shaded area denotes the experimental IVGD resonance energy 13.3 it 0.1 MeV. 
In the right panel the asymmetry energy curves, as functions of the baryon density, are plotted for 
different values of the volume asymmetry 04. 



FIG. 4: The isovector dipole strength distribution (left panel) and transition densities (right panel) 
in ^''^Pb, calculated with a density-dependent effective interaction with 04 = 35 MeV. The proton, 
neutron and total isovector transition densities correspond to the peak at E = 13.3 MeV. 



FIG. 5: The IVGDR excitation energy of ^'^^Pb (upper left panel), and the parameter po of the 
linear density dependence of the nuclear matter asymmetry energy, as functions of the volume 
asymmetry 04. The shaded area denotes the experimental IVGD resonance energy 13.3 it 0.1 MeV. 
In the right panel the asymmetry energy curves, as functions of the baryon density, are plotted for 
different values of the volume asymmetry 04. 



FIG. 6: Average energy gap between the last occupied and first unoccupied major shells in ^ospb, 
as function of the isoscalar parameter 6 (|4^) . The average gaps of neutron states are denoted by 
dots, and those of proton states by squares. 



FIG. 7: Centroid energies of the isoscalar quadrupole Hartree response (upper panel), and the 
ISGQR peak energies calculated in RRPA (lower panel), for five different density-dependent in- 
teractions characterized by the parameter 6 (^6|). As in the previous examples, the calculation is 
performed for ^''^Pb. The shaded area denotes the empirical ISGQR excitation energy in ^'^^Pb: 



10.9 ±0.3 MeV |36|. 



FIG. 8: The isoscalar quadrupole strength distribution (left panel) and transition densities (right 
panel) in ^'^'^Pb, calculated for 6 = 0.93 (see text for description). The vertical bar denotes the 
Ohuj 2^ discrete state. The proton, neutron, and total isoscalar transition densities correspond to 
the ISGQR peak at E = 11.2 MeV excitation energy. 
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